An Equation of State for Anisotropic Solids under Shock Loading 
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An anisotropic equation of state is proposed for accurate extrapolation of high-pressure shock 
Hugoniot states to other thermodynamics states for shocked single crystals and polycrystalline 
alloys. The proposed equation of state represents mathematical and physical generalization of the 
Mie-Griineisen equation of state for isotropic material and reduces to this equation in the limit of 
isotropy. Using an anisotropic nonlinear continuum framework and generalized decomposition of a 
stress tensor [A. A. Lukyanov, Int. J. Plasticity 24, 140 (2008)], the shock waves propagation along 
arbitrary directions in anisotropic solids of any symmetry can be examined. The non-associated 
strength model includes the distortion effect of the yield surface which can be used to describe the 
anisotropic strength differential effect. A numerical calculation showed that the general pulse shape, 
Hugoniot Elastic Limits (HELs), and Hugoniot stress levels for aluminum alloy 7010-T6 agree with 
the experimental data. The results are presented and discussed, and future studies are outlined. 



INTRODUCTION 

The shock waves and high-strain-rate phenomena are 
involved in many physical phenomena, therefore, we are 
interested in understanding the physical properties of 
solids under these non-trivial conditions. Investigation 
of different anisotropic solids (e.g., single crystals and 
polycrystalline alloys) under shock loading has found sig- 
nificant interest in the research community fxj] — fioj] - Mod- 
ern, high-resolution methods to monitoring the stress and 
particle velocity histories in shock waves and equipment 
have been created 



J; numerous investigations into 
the mechanical properties of different classes of materials 
have been undertaken 



lit, [li 



Si]; and numerous phe- 
nomenological as well as microscopic models have been 
developed f 



9|, [30(. However, in spite of a 
perfectly adequate general understanding, experimental 
methodology, and theory, material models do not agree in 
detail, especially for anisotropic solids. For many years, it 
has been assumed that the response of materials to shock 
loading is isotropic, and only recently has anisotropy in 
the shock response attracted the attention of researchers. 
In the early 2000s, Gray III et al. [gj, Q showed that 
under shock loading conditions (one-dimensional strain 
space), the variation of the Hugoniot elastic limit (HEL) 
or the yield strength of annealed zirconium was consis- 
tent with the quasi-static experimental data. To describe 
the anisotropic solids response under shock loading the 
following several factors must be addressed during the 
numerical modelling: (a) large compressions lead to a 
nonlinear behavior - an equation of state (EOS) is re- 
quired, (b) low yield stresses lead to large plastic de- 
formations - an appropriate strength model is required. 
To address these factors, thermodynamically consistent 
framework for modelling the response of single crystal 
and polycrystalline alloys under shock loading was de- 
veloped. This framework, building on the thermody- 
namic approach of Wallace [l[ and continuum framework 
of Johnson couples nonlinear elasticity with non- 



associated anisotropic plasticity within a thermodynam- 
ically consistent numerical incremental formalism [5(. 



SHORTCOMINGS OF THE ISOTROPIC 
EQUATION OF STATE 

The EOS for isotropic materials typically defines the 
pressure as a function of density p (or specific volume, v) 
and specific internal energy e. Experimental shock Hugo- 
niots have been widely used as reference state data for ex- 
trapolating to other high-pressure and high-temperature 
thermodynamic states. The extrapolation has been done 
by using a very popular form of equation of state that is 
used extensively for isotropic solid continua is the Mie- 
Griineisen EOS: 
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where T (v) is the Griineisen gamma, 70 is the initial 
Griineisen gamma, a is the first order volume correction 
to 70, and p is the hydrostatic pressure. Assuming T(v) 
to be a function of volume only JT]), the Mie-Griineisen 
EOS provides a simple extrapolation of a known set of 
Hugoniots data to other sets of thermodynamic states: 
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where Ph is the Hugoniot pressure, p, = 1 is the rel- 

ative change of volume. The shock compression process 
is characterized by the Rankine-Hugoniot jump relations 
lH, [H, [H, [3(|. The Rankine-Hugoniot equations for 
the shock jump conditions includes five variables: pres- 
sure p, particle velocity u p , shock velocity U s , density p 
(or specific volume, v ) and energy e. Hence, an addi- 
tional equation is needed. In many dynamic experiments, 
Up (the particle velocity directly behind the shock) and 
U s (the velocity at which the shock wave propagates 
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through the medium) are measured. The Hugoniot pres- 
sure and a shock velocity U s are a non-linear function 
of particle velocity u p . Following Steinberg's model [l9j |. 
the shock velocity is: 



U s = c + Siu p + S 2 
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(3) 



where c is the velocity curve intercept (can be approxi- 
mated by the bulk speed of sound). Therefore, the Mie- 
Griineisen equation of state with cubic shock velocity as 
a function of particle velocity defines pressure as 
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Qijj is the first generalization of the Kronecker's delta 
symbol. The constructive definition of the tensor «ij is 
based on the fact that stress tensor pa^- is induced in the 
anisotropic medium by the applied isotropic strain tensor 
ydy , i.e. 



pa t j — ——Cijki$ku 



P = -K c e v 



where p is the pressure, e v is the volumetric strain, Sij 
is the Kronecker's delta symbol (unit tensor), Ciju is 
the elastic stiffness matrix, and Kc is the first general- 
ized bulk modulus. The expressions for onj components 
and Kc are presented at the end of Sect. „ Besides, ev- 
erywhere contraction by repeating indexes is assumed. 
Using ([B]), the following expression for total generalized 
"pressure" p* can be obtained: 
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where E is the internal energy per initial specific volume 
e \ 

E = — , K is the classical bulk modulus, Si, S2, S3 
Pa J 

are the intercept of the U-u p curve. Parameters c, Si, 
S2, S3, 70, a represent material properties which define 
its EOS. Parameters have been defined to cover a large 
number of isotropic materials (l9j . 



AN ANISOTROPIC EQUATION OF STATE 

Before discussing an anisotropic equation of state, the 
generalized decomposition of the stress tensor is summa- 
rized in [Io| . The generalized decomposition framework 
will provide a useful point of construction an anisotropic 
equation of state. 



Generalized decomposition of the stress tensor: a-/3 
decomposition 

The definition of pressure in the case of an anisotropic 
solids should be the result of stating that the "pressure" 
term should only produce a change of scale, i.e. isotropic 
state of strain. This should allow the calculation of the 
direct stress ratios that will produce only a change of 
scale The generalized decomposition of the stress 

tensor era is defined as: 
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where P = p*a-ij is the generalized spherical part of 
the stress tensor, S = Sij is the generalized deviatoric 
stress tensor, p* is the total generalized "pressure" and 



where ||a|| = ctijOiij = -f a| 2 + a 2 3 , and finally, the 
expression for the generalized deviatoric part of the stress 
tensor can be rewritten in the following form: 
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For anisotropic materials, the total generalized "pres- 
sure" p* has been expressed [13] as: 
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where p is the pressure related to the volumetric defor- 
mation ([7]) and Pij is the second generalization of the 
Kronecker's delta symbol. The constructive definition of 
the tensor Pij is based on the fact that stress tensor ppij 
is independent of the stress tensor Cijkieu, i-e. their con- 
traction product is zero for any deviatoric strain tensor 
eki, where Cijki is the elastic stiffness matrix. The fol- 
lowing relation describes the functional definition of the 
second order material tensor P^: 



PijCijki = 3KsSti, 



(11) 



where K$ represents the second generalized bulk modu- 
lus. The solution of equations (jTTJ) in terms of Pij compo- 
nents and expression for K$ are presented below. Equa- 
tions (Tl0|) define the correct generalized "pressure" for 
the elastic regime. Note that the generalized decomposi- 
tion of the stress tensor can be applied for all anisotropic 
solids of any symmetry and represents a mathematically 
consistent generalization of the conventional isotropic 
case. The procedure of construction for the tensor aui 
has been defined in (lOj . The elements of the tensor au 
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where Cjj is the elastic stiffness matrix (written in Voigt 
notation). The elements of the tensor (3ki are 
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where are elements of compliance matrix (written 
in Voigt notation). In the limit of isotropy, the pro- 
posed generalization returns to the traditional classical 
case where tensors ay, @ij equal Sij and equations (I10|) 
take the form: 
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where p = — ^ is the classical hydrostatic pressure. 

Also, parameters Kc and K$ were considered as the first 
and the second generalized bulk moduli. In the limit 
of isotropy they reduce to the well-know expression for 
conventional bulk modulus. 



Thermodynamical framework 

During shock loading, the medium undergoes nonlinear 
behavior; therefore an equation of state (EOS) is required 
to describe the medium's response under shock loading 
conditions. It is convenient for use in numerical codes 
to have an analytical form of the EOS. Such an analytic 
form is at best an approximation to the true relationship. 
Thermodynamical definition of the Griineisen parameter 
T(v, T) is the following: 
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where S is the entropy, T is the temperature, C v is the 
heat capacity at constant specific volume. This definition 
may be generalized to a tensor 1^ for the stress-strain 
variables []]]: 
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where £y is the strain tensor and C e is the heat capac- 
ity at constant strains. Using Maxwell's thermodynamic 
relations, the isotropic Mie- Griineisen EOS can be gen- 
eralized for anisotropic media as follows: 



dry 
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where ry is an anisotropic stress tensor produced by the 
isotropic strain state. Also, this EOS (| 19|) is consistent 
with the Wallace's model [l(, where a full anisotropic 
stress was used. It has been shown (see Sect. J that the 
stress tensor patj is produced by the spherical isotropic 
strain tensor edij , hence (|19[) can be re- written as: 



Tij = a l3 T, r = v ( — 



dp 
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where the pressure p is defined by the constitutive equa- 
tions (fTU)) and the isotropic strain. 



An anisotropic EOS 



The equations p0[) define the correct generalized 
"pressure" for the elastic regime. To provide an appro- 
priate description for general anisotropic materials be- 
havior at high pressures, the pressure p related to the 
volumetric deformations is described by p EOS ([4]), ([20)) . 
Therefore, an appropriate description of generalized hy- 
drostatic "pressure" at high pressures has the following 
form: 



p =p 
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(21) 



which also describes correctly the medium's behavior at 
small volumetric strains. To be consistent with the defi- 
nition of the isotropic bulk speed of sound, the following 
definitions of the first c/ and the second C77 bulk speed 
of sound for anisotropic medium are assumed: 
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where the generalized bulk moduli Kc, K$ are defined 
according to (TIB"]) and (fT5|) respectively. Parameters c £ 
[cii,ci], Si, S2, S3, 70, a represent material properties 
which define its EOS f4|. A description of their numerical 
values for AA7010-T6 is shown in Table H 
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TABLE I. EOS data for AA7010-T6 used in analysis. 
Parameter Description AA7010-T6 



c [m/s] Velocity curve intercept 5154 

51 First slope coefficient 1.4 

52 Second slope coefficient 0.0 

53 Third slope coefficient 0.0 
70 Griineisen gamma 2.0 
a First-order volume correction 0.48 



Po [kg/m 3 ] Initial density 2810 
Kc [GPa] Generalized bulk modulus 74.65 



TABLE II. Material properties for AA7010-T6 used in anal- 



ysis. 

Parameter Description AA7010-T6 

F Anisotropy coefficient 0.6898 

G Anisotropy coefficient 0.2873 

H Anisotropy coefficient 0.6824 

Y [MP a] Yield stress 500.0 

X [MPa -1 ] Pressure dependency factor 0.0 

an Tensor aij (11 direction) 0.9976 

Q22 Tensor oiij (22 direction) 1.0029 

«33 Tensor (33 direction) 0.9994 



A NON-ASSOCIATED ANISOTROPIC 
PLASTICITY MODEL 

The main aspects of a phenomenological strength 
model can be characterized by a yield criterion repre- 
senting a surface that separates the elastic and plastic 
regions of the stress space, a flow potential gradient that 
represents the direction of plastic strain rate, a strain 
hardening rule and that plastic flow is incompressible. 



An anisotropic yield surface 



Following the research of Spitzig and Richmond [3l|, 
and Stoughton and Yoon [32j , the mathematically consis- 
tent yield function of a fully anisotropic material based 
on generalized decomposition of the stress tensor is de- 
veloped : 
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where \f r (§ij^j is described by generalized Hill's yield 
functions: 

* 2 (Sij) = F (a 3 S 2 - a 2 S 3 y + G (a x S 3 ~ + 
+H (a 2 2S n - aiiS , 22 ) + 2NS\ 2 + 2LS% 3 + 2MSf 3 

( 24 ) 

where Sij is the generalized deviatoric stress tensor (Si = 
Su, i = 1,2, 3); otij is the generalized Kronecker's symbol 
[1 01 ] {pi = OLn,i = 1,2,3). The material constants Xi 
F, G, H, N, L, M are specified in terms of selected 
initial yield stresses in uniaxial tension, compression, and 
equibiaxial tension. It is important to note that plasticity 
model ([24]) is naturally independent from the generalized 
hydrostatic pressure and therefore, the following equality 
can be written: 

* = * {ay) , Uij = p*a^ + Sij, p* ^ 0, (25) 



where <xy is the stress tensor. In this paper, a uniax- 
ial strain state (one-dimensional reduced mathematical 
formulation in strain space) and the adiabatic approxi- 
mation assumptions are considered for modelling shock 
waves propagation in anisotropic solids. Therefore, ma- 
terial parameters and their numerical values for AA7010- 
T6 are taken as presented in Table [Til 



An anisotropic plastic potential 

A flow potential gradient that represents the direction 
of plastic strain rate is described by the classical Hill's 
anisotropic plastic potential: 



n(<7y-,7y)= F (<7 2 2 -CT 33 ) + G ((T 33 - On) 
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(26) 

and has a strain-rate dependent constant on including 
hardening modulus for plastic potential; 7^ is the back 
stress for plastic potential, D?- is the plastic strain rate 
tensor. It is assumed that the parameters F, G, H, N, 
L, M can be determined from the anisotropy parameters 
R, P, Q23, Q31 and Q12 as for the Hill's yield function 
[33J. The constitutive equations are integrated using the 
tangent stiffness matrix. The numerical values for plastic 
potential parameters were taken from [34j]. 



SIMULATION OF ANISOTROPIC SHOCK WAVE 
PROPAGATION 

The plane shock-wave technique provides a powerful 
tool for st udy ing material properties at different strain 
rates 0, E3, The characteristics as spall pres- 

sure, shock velocity, particle velocity, Hugoniot elastic 
limit, thickness of the spall section, time to spall, and 
free surface velocity of the spall section can be measured 
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and used for the characterization of material dynamic 
response [H, 0, [24 1. 



Description of Experiment 

The flyer plates were launched using a 75 mm bore, 1 
m long and 50 mm bore, and 1 m long single stage gas 
gun the Royal Military College of Science [l5j . The shock 
propagation in the target is monitored using manganin 
stress gauges, placed between target plate and PMMA ( 
poly methylmethacrylate) plate within the target assem- 
bly. The results from the stress gauges were converted by 
N. Bourne and J. Millett to in material (Target) values 
crjvf, using the shock impedances of the target At and 
PMMA Ap, via the well-known relation 



cm 
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(27) 



where ap is the stress gauges values. The 2.5mm thick 
flyer plates of 6082-T6 (dural) were impacted onto the 
targets (test material AA 7010-T6) over the velocity 
range 234 to 895, inducing stresses in the range 2.7 to 
7.2 GPa. The aluminium alloy 6082-T6 was chosen as 
the flyer due to the close similarity in density and wave 
speeds, so that the impact experiments were near sym- 
metrical 34j. The elastic material properties of 7010-T6 
can be found in 34j. Material properties of plates 6082- 
T6 and PMMA can be found in [lj|. 

Using a finite-difference wave propagation code [35| . 
numerical simulations of plate-impact test were per- 
formed with a 5 mm thick target plate, 2.5 mm thick flyer 
plate and 5 mm PMMA plate. Based on the character- 
istics of this plate impact problem, the plates (numerical 
domains), which are used in the numerical simulation, 
are modelled as ID bars 35]. The mesh resolutions were 
sufficient to allow the resolution of all the relevant elas- 
tic and plastic waves in the target and flyer. The stress 
time histories were recorded in the middle of the target 
plate and at the back of the test specimen (the first FD 
element in the PMMA connected to test plate). 



Shock wave propagation in 7010 T6 anisotropic 
aluminum alloy 

The experimental data for AA7010 T6 presented here 
correspond to the plate impact test performed by N. 
Bourne and J. Millett at Royal Military College of Sci- 
ence (published in [34]]) with impact velocities of 450 m/s 
and 895 m/s. Figures Q] and [2] show the comparison be- 
tween experimental data and the numerical simulation 
resulting from the new anisotropic equation of state and 
non-associated anisotropic plasticity model for the longi- 
tudinal and transverse cases. 
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FIG. 1. Back-surface gauge stress traces from plate-impact 
experiments versus numerical simulation of stress (PMMA) 
waves for plate impact test (impact velocity 450 m/s and 895 
m/s) - aluminum alloy AA7010-T6 (Longitudinal Direction). 
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FIG. 2. Back-surface gauge stress traces from plate-impact 
experiments versus numerical simulation of stress (PMMA) 
waves for plate impact test (impact velocity 450 m/s and 895 
m/s) - aluminum alloy AA7010-T6 (Transverse Direction). 



The experimental values, 0.39 GPa and 0.33 GPa for 
elastic response from the longitudinal and short trans- 
verse directions, respectively figures [T] and [21 are in good 
correlation with the modelled values of the HEL longi- 
tudinal - 0.395 GPa and short transverse - 0.333 GPa. 
The errors with respect to the experimental values are 
1.4 % and 0.9 %, respectively to the longitudinal and 
short transverse directions. Further comparison shows 
that the stress pulse width (approximately, 0.81 /is ) and 
release trace are in good agreement with the experiment 
as well (see fig. [T]and[2]). Besides, another important 
characteristic, the arrival time to the HEL and the plastic 
wave velocity are in good correlation with experimental 
data. The main conclusion obtained from these results 
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is that the non-associated anisotropic plasticity model, 
as it stands, is suitable for simulating elastoplastic shock 
wave propagation in anisotropic solids. Different HELs 
are obtained when the material is impacted in different 
directions; their excellent agreement with the experiment 
demonstrates adequateness of the proposed anisotropic 
plasticity model. However, further work is required both 
in the experimental and constitutive modelling areas to 
find a better description of anisotropic material behavior 
and, particulary, Hugoniot stress level. 

CONCLUSIONS 

The EOS proposed in this paper represents a phys- 
ical generalization of the classical Mie-Gruneisen equa- 
tion of state for isotropic materials. Based on the a — [3 
generalized decomposition of stress tensor, the modi- 
fied Mic-Grimciscn equation of state combined with non- 
associated plasticity model forms a system of constitu- 
tive equations suitable for shock wave propagation in sin- 
gle crystals and polycrystalline alloys. The behavior of 
the aluminum alloy 7010-T6 under shock loading condi- 
tions was investigated. A comparison of the experimen- 
tal HELs with numerical simulation shows an excellent 
agreement (maximum error is 1.4%). Furthermore, the 
good agreement of the general pulse shape and Hugoniot 
stress level suggests that the EOS is performing satisfac- 
torily. The significance of the proposed strength model 
includes also the distortion of the yield function shape in 
tension, compression and in different principal directions 
of anisotropy (e.g., 0°, 90°) which can be used to describe 
the anisotropic strength differential effect (anisotropic 
SDE). However, further development of the constitutive 
equations taking into account strain rate sensitivity is 
required. 
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